NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


THESIS 


SWELL  TRANSFORMATION  ACROSS  THE 

CONTINENTAL  SHELF 

by 

Cynthia  Viemes  Tinder 

June  2000 

Thesis  Advisor: 

Thomas  H.C.  Herbers 

Second  Reader: 

Edward  B. Thornton 

Approved  for  public  release;  distribution  is  unlimited. 

[©210  'IJ'TIHTS’  DJZFSamD  4 


20000818  062 


|  REPORT  DOCUMENTATION  PAGE 

Form  Approved  OMB  No.  0704-0188  j 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for 
reviewing  instruction,  searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the 
collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  headquarters  Services,  Directorate  for  Information  Operations  and 
Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget, 
Paperwork  Reduction  Project  (0704-0188)  Washington  DC  20503. 

1.  AGENCY  USE  ONLY  (Leave  blank) 

2.  REPORT  DATE 

June  2000 

3.  REPORT  TYPE  AND  DATES  COVERED 

Master’s  Thesis 

TITLE  AND  SUBTITLE  : 

Swell  Transformation  Across  the  Continental  Shelf 

5.  FUNDING  NUMBERS 

6.  AUTHOR(S) 

Cynthia  Viemes  Tinder 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School 

Monterey,  CA  93943-5000 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESSEES) 

N/A 

10.  SPONSORING  /  MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES  The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official  policy  or 
position  of  the  Department  of  Defense  or  the  U.S.  Government. 

12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited 

12b.  DISTRIBUTION  CODE 

ABSTRACT  (maximum  200  words) 

A  preliminary  analysis  of  swell  decay  across  the  continental  shelf  was  conducted  using  data  from  the  SHOaling  Waves 
Experiment  (SHOWEX).  Six  Datawell  Directional  Waverider  buoys  were  deployed  along  a  cross-shelf  transect,  offshore  of  Duck, 
North  Carolina  during  an  active  hurricane  season  in  the  fall  of  1999.  Estimates  of  the  frequency  spectrum,  E(  f ),  mean  propagation 
direction,  0mean(f\  and  directional  spread,  ae(/),  were  obtained  from  the  auto-,  co-  and  quadrature  spectra  of  the  horizontal  and 
vertical  orbital  displacements  measured  by  the  buoys.  Twelve  cases  were  analyzed  including  eight  cases  dominated  by  remotely 
generated  long  period  swells  arriving  from  various  directions.  The  last  four  cases  were  characterized  by  active  wave  generation  by 
strong  winds  on  the  continental  shelf,  including  the  passages  of  Hurricanes  Floyd  and  Irene.  Comparisons  of  significant  wave 
heights  of  swell  observed  at  all  buoys  show  a  consistent  strong  decay  across  the  shelf  supporting  the  major  role  of  bottom  friction  in 
swell  transformation.  A  decline  in  wave  heights  near  the  coast  was  also  observed  during  the  passage  of  hurricanes,  due  to  limited 
fetch  and  enhanced  bottom  friction  in  shallow  water.  Estimates  of  0mean  show  the  expected  refraction  of  waves  to  shore-normal 
directions  on  the  inner  shelf.  Estimates  of  <j9  generally  do  not  vary  much  across  the  shelf.  However,  a  sharp  increase  of  ae  across  the 
inner  shelf  was  observed  for  directionally  narrow  swells  from  Hurricane  Gert,  possibly  due  to  wave  scattering  from  bottom 
irregularities.  In  all  cases,  CTe  was  at  a  minimum  at  the  spectral  peak  frequency. 

14.  SUBJECT  TERMS  Swell  Transformation,  Continental  Shelf,  Directional  Wave  Spectra,  Wave 
Dissipation,  Wave  Buoy,  Refraction 

15.  NUMBER  OF  PAGES 

72 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

Unclassified 

18.  SECURITY 

CLASSIFICATION  OF  THIS 
PAGE 

Unclassified 

19.  SECURITY 
CLASSIFICATION 
OF  ABSTRACT 
Unclassified 

20.  LIMITATION  OF 
ABSTRACT 

UL 

NSN  7540-01-280-5500  Standard  Form  298  (Rev.  2-89) 


Prescribed  by  ANSI  Std.  239-18 


1 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


u 


Approved  for  public  release;  distribution  is  unlimited 


SWELL  TRANSFORMATION  ACROSS  THE  CONTINENTAL  SHELF 


Cynthia  Viemes  Tinder 
Lieutenant,  United  States  Navy 
B.S.,  United  States  Naval  Academy,  1993 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  METEOROLOGY  AND  PHYSICAL 

OCEANOGRAPHY 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
June  2000 


Author: 


Cynthia  Viemes  Tinder 


Approved  by: 


Thomas  H.C.  Herbers,  Thesis  Advisor 


Edward  B.  Thornton,  Second  Reader 


Roland  W.  Garwood,  Jr.,  Chairman 
Department  of  Oceanography 


m 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


IV 


ABSTRACT 


A  preliminary  analysis  of  swell  decay  across  the  continental  Shelf  was  conducted 
using  data  from  the  SHOaling  Waves  Experiment  (SHOWEX).  Six  Datawell  Directional 
Waverider  buoys  were  deployed  along  a  cross-shelf  transect,  offshore  of  Duck,  North 
Carolina  during  an  active  hurricane  season  in  the  fall  of  1999.  Estimates  of  the  frequency 
spectrum,  E(f),  mean  propagation  direction,  0mecm(f),  and  directional  spread,  ae(/), 
were  obtained  from  the  auto-,  co-  and  quadrature  spectra  of  the  horizontal  and  vertical 
orbital  displacements  measured  by  the  buoys.  Twelve  cases  were  analyzed  including 
eight  cases  dominated  by  remotely  generated  long  period  swells  arriving  from  various 
directions.  The  last  four  cases  were  characterized  by  active  wave  generation  by  strong 
winds  on  the  continental  shelf,  including  the  passages  of  Hurricanes  Floyd  and  Irene. 
Comparisons  of  significant  wave  heights  of  swell  observed  at  all  buoys  show  a  consistent 
strong  decay  across  the  shelf  supporting  the  major  role  of  bottom  friction  in  swell 
transformation.  A  decline  in  wave  heights  near  the  coast  was  also  observed  during  the 
passage  of  hurricanes,  due  to  limited  fetch  and  enhanced  bottom  friction  in  shallow  water. 
Estimates  of  6mean  show  the  expected  refraction  of  waves  to  shore-normal  directions  on 
the  inner  shelf.  Estimates  of  Oe  generally  do  not  vary  much  across  the  shelf.  However,  a 
sharp  increase  of  Cq  across  the  inner  shelf  was  observed  for  directionally  narrow  swells 
from  Hurricane  Gert,  possibly  due  to  wave  scattering  from  bottom  irregularities.  In  all 
cases,  Ce  was  at  a  minimum  at  the  spectral  peak  frequency. 
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I. 


INTRODUCTION 


Swell  characteristics  transform  dramatically  between  the  open  ocean  and  the 
beach.  When  the  water  depth  becomes  comparable  to  the  swell  wavelength,  the  bottom 
starts  to  influence  swell  propagation  by  causing  changes  in  the  wave  phase  and  group 
speeds. 

As  swells  begin  to  shoal  the  group  speed  initially  increases,  and  it  follows  from 
the  conservation  of  energy  that  the  wave  amplitude  decreases  by  as  much  as  10%  in 
intermediate  water  depths  (kh  ~  1,  where  k  is  the  wave  number  and  h  is  the  water  depth). 
In  shallower  water,  the  group  speed  decreases  with  decreasing  depth,  causing  an  increase 
in  wave  amplitudes  near  the  coast. 

The  phase  speed  changes  associated  with  depth  variations  cause  refraction 
analogous  to  Snell’s  law  in  geometric  optics.  Monochromatic  wave  trains  propagating 
over  a  shelf  with  straight  and  parallel  depth  contours,  refract  Such  that  the  crests  become 
more  parallel  to  the  beach.  In  situations  with  complex  bathymetry,  refraction  can  cause 
large  wave  height  variations.  Wave  propagation  directions  generally  refract  towards 
shallower  regions,  typically  causing  a  convergence  region  of  larger  waves  near  shoals, 
headlands  and  ridges,  and  divergence  regions  of  smaller  waves  in  bays  and  near 
submarine  canyons  (e.g.,  Munk  and  Traylor,  1947).  Significant  refraction  effects  occur 
as  waves  travel  over  a  large  submerged  topographic  feature,  but  also  when  waves  move 
through  a  large  shallow  coastal  region  with  smooth  undulating  topography  (e.g.,  Graber 
et  al.,  1990;  O’Reilly  and  Guza,  1993). 
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Shoaling  and  refraction  are  linear  processes  that  are  well  understood  theoretically. 
By  applying  an  energy  balance  along  ray  trajectories,  the  evolution  of  wave  spectra  can 
be  evaluated  for  arbitrary  two-dimensional  continental  shelf  topography  (e.g.,  Munk  and 
Arthur,  1951;  Longuet-Higgins,  1957;  Le  Mehaute  and  Wang,  1982). 

Although  swell  can  cross  entire  ocean  basins  with  minor  energy  losses  (Snodgrass 
et  al.,  1966),  significant  dissipation  of  swell  energy  takes  place  on  a  shallow  continental 
shelf.  There  are  several  mechanisms  through  which  the  bottom  can  remove  energy  from 
the  waves  and  reduce  wave  heights.  Shemdin  et  al.  (1980)  evaluated  the  relative 
importance  of  bottom  friction  (usually  most  dominant  in  sandy  bottoms  with  0.1 -0.4  mm 
grain  diameter),  percolation  through  the  bottom  (coarser  sediments  with  diameter  >  0.5 
mm),  and  wave-induced  bottom  motion  (fine  sediment).  The  sediment  type,  bottom 
microtopography,  in  particular,  the  presence  of  ripples,  and  local  currents  play  important 
roles  in  the  damping  rates,  all  of  which  are  often  unknown. 

On  the  sandy  continental  shelf  considered  in  this  study,  bottom  friction  is  believed 
to  be  the  dominant  swell  decay  mechanism.  Bottom  friction  is  dependent  on  both  sand 
ripple  characteristics  and  the  prevailing  sea-state  conditions  at  the  site.  Tolman  (1994) 
reviewed  the  important  feed  back  between  ripple  formation  and  wave  energy  dissipation. 
Orbital  sand  ripples,  initiated  when  the  orbital  wave  motion  near  the  bottom  is  strong 
enough  to  move  sediment,  greatly  enhance  the  bottom  roughness  (Grant  and  Madsen, 
1982).  Tolman’s  (1994)  model  of  wave-generated  bedforms  and  associated  damping  of 
swell  across  a  continental  shelf  suggest  that  moderately  energetic  swell  shoaling  on  a 
continental  shelf  can  generate  sand  ripples  that  cause  strong  damping  of  the  swell. 
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Few  measurements  of  the  attenuation  of  swell  in  shallow  water  have  been 


reported.  The  JONS  WAP  experiment  examined  the  swell  decay  observed  over  a  160- 
km-long  cross-shore  transect  of  various  instruments  deployed  off  the  coast  of  Germany 
(Hasselmann  et  al.,  1973).  The  observed  damping  rates  did  not  agree  with  generally 
accepted  formulations  of  bottom  friction  and  suggested  that  other  attenuation 
mechanisms  such  as  back-scattering  of  swell  from  small  scale  bottom  irregularities  may 
be  important  (Long,  1973).  Young  and  Gorman  (1995)  conducted  a  2-week-long 
experiment  using  a  cross-shore  transect  of  various  types  of  instruments  spanning  the 
continental  shelf  on  the  southern  coast  of  Australia.  These  measurements  were  used  in 
conjunction  with  the  spectral  wave  model  WAM  (The  WAMDI  Group,  1988)  to  estimate 
the  contribution  of  bottom  friction  to  the  observed  overall  decrease  in  swell  energy  across 
the  shelf.  The  results  suggested  that  the  strong  attenuation  of  energetic  southern  ocean 
swells  observed  at  shallow  instrument  sites  is  primarily  caused  by  bottom  friction  with 
estimated  drag  coefficients  that  are  much  larger  than  the  values  used  in  wave  prediction 
models.  Forristall  and  Reece  (1985)  measured  wave  attenuation  and  bottom  processes 
near  the  Mississippi  Delta  in  the  Sea  Wave  Attenuation  Measurement  Program 
(SWAMP).  They  suspected  that  the  soft  bottom  near  river  deltas  causes  strong  swell 
attenuation.  The  observed  attenuation  was  a  strong  function  of  the  deep-water  significant 
wave  height  and  a  weak  function  of  the  wave  frequency.  Refraction  and  shoaling 
accounted  for  the  changes  in  the  directional  spectra  when  the  wave  height  was  small. 

With  larger  wave  heights  (>  3  m),  significant  attenuation  was  observed,  in  particular  for 
frequencies  <0.10  Hz. 
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Herbers  et  al.  (2000)  studied  the  transformation  of  swell  across  the  continental 
shelf  off  the  coast  of  North  Carolina.  The  results  suggested  that  on  the  outer  shelf  energy 
was  lost  only  by  the  lower  frequencies  of  the  spectrum,  possibly  because  high-frequency 
components  with  relatively  shorter  wavelengths  were  unattenuated  at  the  bottom.  At  the 
shallower  sites,  observed  spectral  levels  were  lower  than  predicted  spectral  levels  at  all 
frequencies,  suggesting  significant  damping  of  all  components  of  the  spectrum.  The 
apparent  damping  is  strongest  near  the  peak  frequency  causing  a  slight  broadening  of  the 
observed  spectrum  at  the  closest  nearshore  site.  The  inferred  friction  factor  agreed 
roughly  with  the  results  from  Young  and  Gorman’s  study,  and  was  much  larger  than 
those  commonly  used  in  wave  prediction  models  (e.g.,  The  WAMDI  Group,  1988). 

New  observations  of  the  evolution  of  swell  across  a  continental  shelf  are 
presented  here.  Six  Datawell  Directional  Waverider  buoys  were  deployed  along  a  100 
km  cross-shelf  transect  off  the  coast  of  North  Carolina  for  3  months  in  1999  as  part  of  the 
SHOWEX  (SHoaling  Waves  Experiment).  In  contrast  to  previous  experiments  that 
provided  estimates  of  swell  height  decay  across  the  shelf,  the  detailed  directional  wave 
measurements  collected  in  SHOWEX  also  allow  for  an  analysis  of  the  transformation  of 
directional  wave  properties  across  the  shelf.  An  unusually  active  hurricane  season  caused 
by  La  Nina  conditions  provided  numerous  case  studies  of  energetic  long  period  swell. 

The  field  experiment  is  described  in  Chapter  II,  followed  by  a  description  of  the  data 
analysis  procedures  in  Chapter  III.  Twelve  case  studies  are  presented  in  Chapter  IV, 
followed  by  an  analysis  of  observed  cross-shelf  variations  in  Chapter  V.  Chapter  VI 
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discusses  the  frequency-directional  spectra,  and  the  results  are  summarized  in  Chapter 
VII. 
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II.  FIELD  EXPERIMENT 


The  field  data  used  in  this  study  were  collected  as  part  of  the  ONR  SHOaling 
Waves  Experiment  (SHOWEX)  conducted  offshore  of  the  U.S.  Army  Corps  of 
Engineers  Field  Research  Facility  (FRF)  near  Duck,  North  Carolina,  from  mid-September 
to  mid-December,  1999.  The  site  is  located  about  midway  between  Cape  Hatteras  and 
Chesapeake  Bay,  along  a  coast  consisting  of  a  series  of  relatively  straight  barrier  islands 
with  sandy  beaches  that  are  fully  exposed  to  the  Atlantic  Ocean. 

The  continental  shelf  is  oriented  approximately  north-south,  and  widens  from 
Cape  Hatteras  toward  Chesapeake  Bay.  The  continental  shelf  is  about  100  km  wide  and 
only  20  -  50  m  deep  across  most  of  the  study  site  (Figure  1).  Bathymetry  is  characterized 
by  a  gently  sloping  inner  shelf  (about  1 :250),  extending  into  a  wide  mid-shelf  region  with 
an  even  lower  slope  of  about  1 :2500  and  large  (0(5  km))  ridge-like  features  that  have 
vertical  scales  of  0(5  m).  The  continental  shelf  widens  from  Cape  Hatteras,  North 
Carolina,  and  curves  to  the  northeast  toward  Nova  Scotia  and  the  Grand  Banks. 

Six  Datawell  Directional  Waverider  buoys  (XI  to  X6)  were  deployed  along  a 
cross-shelf  transect  extending  almost  due  east  from  the  Duck  beach  to  the  shelf  break 
(Figure  2).  The  buoys  were  deployed  in  water  depths  of  21m  (XI,  located  5.35  km  from 
shore);  24  m  (X2, 12.78  km);  26  m  (X3, 22.17  km  from  shore);  33  m  (X4, 37.74  km  from 
shore);  39  m  (X5,  58.22  km  from  shore);  and  193  m  (X6, 86.57  km  from  shore)  (Table 
!)• 
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Buoy 

Distance 
From  Shore 

Depth 

Latitude 

Longitude 

Dates 

Deployed 

XI 

5.4  km 

21  m 

36°  13.62’ N 

75°  42.28’ W 

9/13/99- 

12/13/99 

X2 

12.8  km 

24  m 

36°  13.62’ N 

75°  36.77’  W 

9/1 1/99  - 
12/12/99 

X3 

36°  12.27’  N 

75°  29.91 ’W 

9/10/99- 

12/13/99 

X4 

9/10/99  - 
12/13/99 

X5 

58.2  km 

39  m 

36°  07.79’  N 

9/13/99  - 
10/24/99 
10/28/99- 
12/13/99 

X6 

86.62  km 

197  m 

36°  04.98’  N 

75°  47.47’  W 

9/13/99- 

9/21/99 

86.57  km 

193  m 

36°  04.97’  N 

75°  47.50’ W 

10/15/99- 

12/5/99 

Table  1 .  Locations  and  time  periods  of  Datawell  1 
(XI  to  X6)  deployments 

Directional  Waverider  buoys 

Buoy  X5  broke  loose  on  October  24,  and  was  redeployed  on  October  28.  Buoy  X6  broke 
loose  on  September  21,  was  redeployed  on  October  15,  and  lost  again  from  December  5 
through  the  end  of  the  experiment.  All  other  buoys  were  operational  throughput  the 
experiment. 

The  Datawell  Directional  Waverider  buoys  are  surface  following  directional  wave 
buoys  (about  1  m  diameter)  which  measure  horizontal  and  vertical  orbital  displacements 
at  the  sea  surface  from  which  wave  height  spectra  and  directional  wave  properties  are 
determined.  They  contain  a  heave-pitch-roll  sensor  (Hippy-40),  three  axis  fluxgate 
compass,  two  fixed  ‘x’  and  *y’  accelerometers,  a  temperature  sensor  and  a 
microprocessor.  The  direction  measurement  is  based  on  the  translation  principle,  using 
the  measured  horizontal  accelerations  rather  than  surface  slopes  used  in  earlier  pitch-and- 
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roll  type  buoys.  Therefore,  the  measurement  is  independent  of  buoy  roll  motions  and 
thus,  a  small  spherical  buoy  can  be  used.  From  the  accelerations  and  tilts  measured  in  the 
moving  ‘buoy  reference  frame,’  the  accelerations  along  the  fixed  horizontal,  north  and 
west  axes  are  calculated.  All  three  accelerations  (vertical,  north  and  west)  are  digitally 
integrated  to  get  filtered  displacements  with  a  high  frequency  cut-off  at  0.6  Hz.  To  save 
transmitting  power,  the  real  time  data  are  compressed  to  vertical,  north  and  west 
displacements.  Every  half  hour,  Fast  Fourier  Transforms  of  8  series  of  256  data  points 
(200  sec)  are  used  to  form  spectra  and  cross-spectra  with  16  degrees  of  freedom. 
Frequency  resolutions  of  0.005  Hz  in  the  range  0.025  to  0.1  Hz,  and  0.01  Hz  in  the  range 
0.1  to  0.58  Hz  are  used.  For  each  frequency  bandwidth,  the  on-board  microprocessor 
computes  energy  density,  mean  direction,  directional  spread,  and  the  normalized  second 
harmonic  of  the  directional  distribution  (discussed  in  Chapter  III).  The  Directional 
Waverider  buoy  transmits  the  following  information:  vertical,  north,  and  west 
displacements,  spectral  density,  processed  directional  parameters,  sea  surface  temperature 
(SST)  and  other  system  status  with  parity  bits  for  error  checking  purposes.  The  buoy 
data,  framed  in  64  bit  vectors,  is  transmitted  via  HF  in  the  27-40  MHz  band  continuously. 
The  processed  spectral  and  directional  data  is  transmitted  8  times  to  allow  for  reliable 
data  return  even  in  marginal  transmission  conditions. 

A  single  point  vertical  mooring  (Figure  3)  ensures  sufficient  symmetrical 
horizontal  buoy  response  even  for  small  motions  at  low  frequencies.  The  Directional 
Waverider  is  fitted  with  5  kg  chain  ballast  to  provide  stability  when  only  a  vertical 
mooring  force  is  present  (free  floating  or  moored  in  shallow  water).  The  low  stiffness  of 
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the  30  m  rubber  cord  allows  the  buoy  to  follow  waves  up  to  40  m  in  height.  High  shock 
accelerations  are  absorbed  by  the  fluid  surrounding  the  accelerometer  and  the  shockproof 
mounting  of  the  accelerometer  housing,  electronics  and  batteries. 

The  four  buoys  located  closest  to  shore  (XI  to  X4)  transmitted  their  respective 
data  via  an  HF  radio  link  to  receivers  mounted  on  a  50  m  high  tower  at  the  FRF,  while 
the  outer  two  buoys  (X5  and  X6)  that  were  out  of  HF  range,  used  internal  data  loggers. 

Significant  weather  associated  with  La  Nina  conditions  played  an  important  role 
in  the  formation  of  long  period  energetic  swell  that  propagated  across  the  continental 
shelf  during  the  fall  of  1 999.  Hurricane  season  peaked  in  September  and  October  with 
Hurricanes  Floyd,  Gert,  Irene  and  Jose  (Figure  4).  On  September  17,  the  eye  of 
Hurricane  Floyd  passed  within  100  km  of  the  buoy  transect  generating  large  seas  on  the 
continental  shelf  and  causing  massive  flooding  and  destruction  in  eastern  North  Carolina. 
Hurricanes  Gert  and  Jose  remained  offshore  but  radiated  long  period  swell  across  the 
transect  of  buoys.  Hurricane  Irene  made  landfall  across  western  Cuba  and  southeastern 
Florida  and  continued  northeastward,  following  the  Gulf  Stream  and  remained  offshore 
of  the  Carolinas.  Further  discussion  on  these  storms  is  found  in  the  case  studies  in 
Chapter  IV. 
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III. 


DATA  ANALYSIS 


The  directional  properties  of  ocean  surface  waves  are  often  described  by  the  mean 
wave  direction,  8mem  (/) ,  and  the  directional  spread,  <J6{f) ,  as  functions  of  frequency. 
These  variables  are  closely  related  to  the  directional  distribution  of  wave  energy, 

S(9;  f) :  0mean  (/)  is  the  mean  of  this  distribution,  defining  a  dominant  wave  propagation 
direction  at  each  frequency,  whereas  cre(f)  represents  the  standard  deviation  of  the 
distribution  and  characterizes  the  half-width  of  the  directional  distribution.  Estimates  of 
9mean  (/)  and  <Jd  (/)  are  readily  obtained  from  heave-pitch-and-roll  or  translational 
directional  wave  buoys  (e.g.,  Longuet-Higgins  et  al.,  1963;  Mitsuyasu  et  al.,  1975;  Long, 
1980,  Kuik  et  al.,  1988).  Longuet-Higgins  et  al.  (1963)  showed  that  the  first  four 
directional  Fourier  coefficients  of  S(9;  f)  can  be  derived  from  the  auto-,  co-  and 
quadrature  spectra  of  the  three  basic  heave  and  slope  (or  horizontal  displacement)  signals 
of  the  buoy.  For  a  Datawell  Directional  Waverider  buoy  that  measures  the  horizontal 
(x,y)  and  vertical  (z)  displacements  of  water  particles  at  the  surface,  the  Fourier  moments 
can  be  expressed  as: 


<*(/)•  cos  (0)  £«?;/)  = 


QJJ± 

{C=(/)[C„(/)  +  Cw(/)]},/2 


(1) 
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(2) 


W)«  jd0sm(0)S(0;f) 


QM} 

{cc(/)[c„(/)  +  c„(/)]},/2 


<h<J)m  jd&  cos  (26)  S (0;  f)  = 


Q(/)-cff(/) 


(3) 


&*(/)■  | d9  sin  (20)  S(0;f)  = 


2CM) 

CJJ)+C„<J) 


(4) 


where  Cy  and  Qy  denote  the  co-  and  quadrature  spectra  of  components  i  and  j 
respectively. 

For  narrow  S(0) ,  0mean  and  oe  can  be  approximated  by  the  first-order  moments, 
ax  and  bx  (e.g.,  Kuik  et  al.,  1988): 


jd9an(9-9mm)S(9)  =  b,  cos (»_)-«,  sin(9_)=0  (5) 


of  -  J<W(2[1-  cos(0-<U,)])S(0)=2[l-(a,  cos(0+*l  sin  (*_))].  (6) 


For  small  9-6 ^ 
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sin  (0-0_)  -  (0-3«»)[l-i(0-*U)2] 

2[i-cos 

and  thus,  0mam  and  a6  represent  a  mean  value  and  standard  deviation  of  the  directional 
distribution  of  wave  energy  with  an  O  (0-0mean  )2  bias. 

A  disadvantage  of  this  formulation  is  the  sensitivity  of  cre  to  directionally 
opposing  wave  components.  A  small  amount  of  wave  energy  traveling  in  a  direction 
opposing  the  dominant  swell  system  (e.g.,  reflection  of  waves  from  shore  or  a  weak 
secondary  swell  peak)  causes  a  large  increase  in  thecr^  estimate  that  is  not  related  to  the 
beam  width  of  the  dominant  swell  .(Herbers  et  al,  1999). 

Alternatively,  0mean  and  g6  can  be  defined  in  terms  of  the  second-order  moments 
a2  and  b2  (e.g.,  Herbers  et  al.,  1999): 


)de  isinP^-fU)]^)"  \[b,  cos(2 sin(20meaJ1)]  =  O  (7) 


o>2  -  )d$  |{[l-cos  P(S-tf_)]>  S(S)  =  i[l-(a2  cos  (2 .«_)+*,  sin  <2»_))](8) 


For  small  0-0me3D 


13 


I  sin  [2(^-<9_  )]  -  )[1  - 1  (e-e^?} 

\  0  -cos[2 (6-emm )] } » («-*_  )![1  -  5  (e-emm  )2] 

and  therefore  the  second-order  approximations  of  0mean  and  cre  have  larger  biases  than 
the  first-order  approximations,  Eqs  5  and  6.  However,  Eqs  7  and  8  are  much  less 
sensitive  than  Eqs  5  and  6,  to  a  weak  directionally  opposing  wave  system,  because  the 
double-angle  argument,  2(0-#mean),is  the  same  for  waves  propagating  in  opposite 
directions.  Unlike  the  first-order  estimate  of  0mtm  (Eq  5),  the  second-order  estimate  has  a 
180°  ambiguity. 

In  the  present  analysis,  estimates  of  0mtm  and  crd  were  obtained  from  both  the 
first-  and  second-order  moments.  Because  the  two  methods  have  different  errors, 
intercomparison  of  0mean  and  <rd  estimates  based  on  a,,  6,  and  a2,  b2  (discussed  below) 
provides  a  useful  check  of  the  accuracy  of  the  estimates.  A  magnetic  correction  of 
-10.75°  W  was  applied  to  0mvm  (/),  and  the  resulting  direction  was  transformed  to 
nautical  convention,  that  is  the  transformed  0mean  (/)  is  the  direction  the  waves  come 
from,  measured  counterclockwise  from  true  north.  In  addition  to  0raean  and  <Je  estimates 
as  functions  of  frequency,  bulk  estimates  were  computed  for  the  frequency  range  0.025  to 
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0.25  Hz.  The  estimates  are  based  on  Eqs  5-8  using  energy-weighted  average  values  of 


the  moments,  av  bx  and  a2,  b2: 


025Bz 

J  <>,(/)£(/)  if 

n  0.025  Hz _ 

W1  total  0.25  Hz 

J  EU)if 

0.025  Hz 


0.25  Hz 

J  4 (/)£(/)  if 

u  —  0.025  Hz _ 

0.25  Hz 

/  £(/)4T 

0.025  Hz 


025  Hz 

|  <%(/)£(/)  ¥ 

n  0.025  Hz _ 

U2<a,oi  025  Hz 

J  E(f>¥ 

0.025  Hz 


025  Hz 

j  b2(j)E<_r>df 

L  0.025  Hz _ 

U2total  025  Hz 

J  £io<r 

0.025  Hz 


(9) 


(10) 


(11) 


(12) 


where  E(f)  [  =  C_(/)]  is  the  wave  frequency  spectrum. 

A  preliminary  analysis  of  the  bulk  wave  and  wind  parameters  during  the  entire 
experiment  was  conducted  to  select  specific  events  for  further  analysis.  For  each  case 


covering  time  periods  of  12  or  3  hours,  the  spectra  of  the  half-hour  data  records  were 
ensemble  averaged,  and  directional  moments,  0mtm  (/)  and  cre(/),  were  computed  based 
on  the  ensemble  averages  of  the  moments  (a,,  bx,  a2,  and  b2)  weighted  by  E(f).  The 
resulting  E{f),  0mean  (/),  and  <Jg{f)  have  approximately  384  degrees  of  freedom. 
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IV. 


CASE  STUDIES 


The  times  of  the  selected  twelve  case  studies  are  listed  in  Table  2,  together  with 
the  significant  wave  height: 


1 0.25  Hz 

H,=4  J  E(f)df,  (13) 

V  0.025  Hz 

the  peak  period,  Tp,  (the  wave  period  corresponding  to  the  peak  of  the  spectrum,  E{f )), 
and  the  bulk  Gmesa  and  crg  parameters  (calculated  using  Eqs  5-12)  all  from  the  offshore 
buoy,  and  the  wave  height  decay  across  the  shelf  quantified  by  the  ratio  HS(X1)/HS(X6). 

Eight  of  the  cases  were  collected  during  light  wind  conditions  on  the  shelf  when 
the  wave  field  was  dominated  by  remotely  generated  long  period  (Tp=12-17  s)  swell. 
Since  it  takes  several  hours  for  the  swell  to  traverse  the  shelf,  12-hour  periods  were  used 
in  the  analysis  with  relatively  steady  significant  wave  heights,  peak  periods  and  mean 
propagation  directions. 

Six  of  the  swell  cases  are  “easterly  swell”  events  with  mean  offshore  directions 
within  +/-  45°  of  090°  T.  The  easterly  swell  cases  include  forerunners  from  Hurricane 
Floyd  and  arrivals  from  Hurricane  Gert.  Hurricane  Gert  was  a  category  IV  hurricane  that 
was  tracking  to  the  northwest  and  slowly  turned  north-northeastward  on  September  21 . 
Gert’s  closest  point  of  approach  was  about  1100  km  to  the  southeast  of  the  buoy  transect 
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(Figure  4).  Gert  produced  the  longest  period  (Tp=  15-17  s)  and  directionally  narrowest 
( a6  =20-24°)  swells  of  the  experiment. 


EASTERLY  SWELL 

Hs  at  X6 

Tp  at 
X6 

©mean  at 

X6 

CTe  at 
X6 

09/14  00  EST  to  09/14  12  EST 
(FLOYD) 

2.90  m 

13.33  s 

0.57 

09/19  12  EST  to  09/20  00  EST 
(GERT) 

2.78  m 

15.38  s 

124.3°T 

0.47 

09/20  21  EST  to  09/21  09  EST 
(GERT) 

2.54  m 

16.67  s 

127.0°  T 

0.55 

X5: 
12.50  s 

X5: 

73 .4°  T 

X5: 

42.5° 

(X1)/(X5): 

0.63 

11/18  03  EST  to  11/18  15  EST 

1.19m 

14.29  s 

54.5°  T 

0.66 

11/22  18  EST  to  11/23  06  EST 

2.07  m 

11.76  s 

97.0  6  T 

0.67 

LARGE  OBLIQUE  ANGLE 
SWELL 

10/23  02  EST  to  10/23  14  EST 
(JOSE) 

1.70  m 

12.50  s 

163.4°T 

m 

0.51 

12/01  21  EST  to  12/02  09  EST 

4.60  m 

0.50 

STORM  WAVES 

09/16  09  EST  to  09/16  12  EST 
(FLOYD) 

7.63  m 

12.50  s 

HI 

0.62 

10/18  00  EST  to  10/18  03  EST 
(IRENE) 

4.30  m 

0.62 

11/02  20  EST  to  11/02  23  EST 

5.00  m 

10.00  s 

167.9°T 

31° 

0.42 

11/1 123  EST  to  11/12  02  EST 

3.58  m 

■IWIKtiia 

Hi fiH 

0.77 

Table  2.  Twelve  case  studies  classified  as:  “Easterly  Swell”,  “Large  Oblique  Angle 
Swell”  and  “Storm  Waves”  listed  at  the  offshore  buoy,  X6,  are:  Significant  wave  height, 
Hs;  peak  period,  TP;  bulk  mean  propagation  direction,  0mean;  and  bulk  directional  spread, 
<7e.  Also  listed  is  the  wave  height  decay  between  XI  and  X6. 

In  the  remaining  two  cases,  swell  arrived  at  very  large  angles  to  shore-normal 
resulting  in  strong  refraction  effects.  The  two  “large  oblique  angle  swell”  cases  include 
southerly  swell  from  Hurricane  Jose  and  a  northerly  swell  case.  Hurricane  Jose  remained 
well  offshore  and  followed  a  somewhat  similar  track  to  Hurricane  Gert,  traveling  to  the 
north-northwest  until  passing  the  eastern  tip  of  Puerto  Rico  before  recurving  to  the  north- 
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northeast  (Figure  4).  Jose’s  closest  point  of  approach  was  900  km  of  the  buoy  transect  to 
the  southeast.  The  southeasterly  swell  from  Jose  propagated  over  the  steep  slope  of  the 
continental  shelf  located  to  the  southeast  of  the  buoy  transect  (Figure  1)  near  Cape 
Hatteras  Point.  In  contrast,  waves  arriving  from  northerly  directions  travel  a  long 
distance  over  the  continental  shelf,  which  widens  and  curves  to  the  northeast. 

In  addition  to  the  8  swell  cases,  4  storm  cases  with  strong  winds  and  active  wave 
generation  on  the  continental  shelf,  were  selected  for  analysis.  Because  of  the  more 
rapidly  varying  conditions,  3 -hour  periods  at  the  peak  of  the  storm  events  were  analyzed. 
These  cases  include  Hurricanes  Floyd  and  Irene,  the  most  significant  storm  events  during 
the  entire  3 -month  experiment.  Their  tracks  passed  within  100-200  km  from  the  buoy 
transect  and  produced  large  waves  (maximum  Hs  of  4.3  -  7.6  m)  on  the  continental  shelf. 
Figure  5  displays  the  associated  winds  and  Hs  of  Hurricanes  Floyd  and  Irene.  The  peak 
winds  of  35  m/s  for  Floyd  occurred  around  2000  EST  on  September  16  (well  after 
Hurricane  Floyd  made  land  fall)  while  the  peak  Hs  at  buoy  X6  were  observed  about  2 
hours  later  (Figure  5a).  Since  Floyd  passed  to  the  west  of  the  buoy  transect  (Figure  4), 
the  strongest  winds  were  observed  closest  to  shore  as  indicated  by  the  difference  in 
maximum  wind  speeds  between  FRF  and  NOAA  Buoy  44014  observations.  The  NOAA 
Buoy  44014  was  located  about  30  km  northeast  of  the  buoy  transect  (Figure  1). 

Hurricane  Irene  made  landfall  over  southeastern  Florida,  but  remained  offshore  of  the 
Carolinas,  and  passed  within  250  km  of  the  buoy  transect  (Figure  4).  As  Irene  passed  to 
the  southeast  of  the  transect,  the  dominant  waves  continued  to  arrive  from  the  south- 
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southeast  even  though  the  strong  local  winds  (on  the  west  side  of  the  hurricane)  were 
from  the  northwest  (Figure  5b). 

The  other  two  storm  events  selected  were  local  storms  (a  cold  frontal  system  on 
November  2  and  a  nor’easter  on  November  1 1-12)  that  produced  comparably  large  waves 
on  the  shelf.  In  the  Floyd,  Irene  and  the  November  2  storm  events,  when  the  dominant 
seas  observed  at  X6  arrived  from  the  south,  strong  refraction  and  fetch-limited  conditions 
on  the  inner  shelf  were  expected.  The  nor’easter  case  shows  very  different  conditions  of 
fetch-unlimited  seas  generated  by  northeasterly  winds. 

The  November  2  storm  moved  slowly  across  the  study  site  with  wind  speed 
observations  from  the  FRF  at  a  maximum  of  18  m/s  at  1200  EST  while  the  NOAA 
buoy’s  peak  wind  speed  occurred  about  10  hours  later,  during  the  3-hour  period  selected 
(Figure  6a).  In  contrast,  the  nor’easter  storm  quickly  passed  through  the  study  area  with 
the  peak  wind  observations  at  the  NOAA  buoy  occurring  3-4  hours  later  than  at  the  FRF. 
Hs  was  at  a  maximum  at  buoy  X6  about  2  hours  after  the  peak  wind  observations  from 
the  FRF  station  (Figure  6b). 

For  each  case  the  transfer  function,  G(f),  was  examined  between  the  horizontal 
and  vertical  displacements,  to  verify  the  buoy  response: 


W)  =  , 


\CM)+C^J) 


cB(f) 


(14) 


For  linear  waves  G(f)  =  1  /  tanh(  k(f)h),  where  h  is  the  water  depth  and  the 
wavenumber,  k,  is  given  by  the  dispersion  relation  (2n  f  )2  =  gk  tanh(&/2).  Hence,  the 
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normalized  transfer  function  G'(f)  [  =  tanh  (k(f)  h )  G(f)]  is  equal  to  1  for  linear 
,  waves.  Nonlinearities  in  both  the  wave  field  and  buoy  response  will  result  in  deviations 
of  G  ’  from  1 .  Bulk  values  of  G  and  G  ’  integrated  over  the  spectra  are  expressed  as: 


G  = 


025  Hz 

f  CJf)+C„(f)df 


t 


0.025  Hz _ • 

0.25  Hz 

1  c„(f)df 

0.025  Hz 


(15) 


0.25  Hz 

I  taWKf)  *)  [C„  (/)  +  C„  ( f)W 

0.025  Hz _ 

0.25  Hz 

1  CJS)if 

0.025  Hz 


(16) 


Figure  7a  shows  the  observed  ratio,  G,  for  each  of  the  twelve  cases.  The  observed  ratios 
are  close  to  1  at  the  deep  water  (buoy  X6),  but,  for  all  cases,  G  generally  increases  from 
the  outer  shelf  to  the  inner  shelf.  This  trend  is  consistent  with  the  transformation  of  the 
wave  orbital  motion  from  deep  (circular  orbits,  G  —  1)  to  shallow  (ellipsoidal  orbits,  G  > 
1)  water.  The  normalized  ratios  vary  between  1  to  1.2  (Figure  7b),  showing  a  slight  but 
systematic  positive  bias  that  increases  with  distance  from  shore,  possibly  due  to  nonlinear 
effects. 

The  accuracy  of  the  directional  moments,  0mean  and<re ,  is  verified  in  Figure  8  by 
comparing  the  estimates  based  on  the  second-order  approximations  (a2,  b2)  for  each  of 
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the  twelve  cases  listed  in  Table  2,  with  their  respective  estimates  based  on  first-order 
approximations  (ax,  bx ).  The  bulk  0maa  values  obtained  with  the  two  different 
approximations  are  in  excellent  agreement  (Figure  8a).  However,  the  ae  values  are  well 
correlated  but  estimates  based  on  (a2,  b2)  values  are  generally  about  15-25%  smaller 
than  estimates  based  on  (a,,  bx )  (Figure  8b).  These  differences  are  qualitatively 
consistent  with  the  larger  negative  bias  expected  for  <je  estimates  based  on  the  second- 
order  moments  discussed  earlier.  The  estimates  of  0mt3B  andcre  used  in  Chapters  V  and 
VI  for  the  preliminary  analysis  of  the  twelve  cases  are  based  on  the  less-biased  first-order 
approximations  (a,,  bx ). 
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V.  CROSS-SHELF  VARIATIONS 


Variations  in  wave  properties  across  the  continental  shelf  can  result  from  different 
attenuation  mechanisms  discussed  in  Chapter  1 .  Bottom  friction,  refraction  of  waves 
arriving  at  oblique  angles,  and  fetch  limitations  of  wave  growth  on  the  continental  shelf 
all  contribute  to  a  decrease  in  wave  height  across  the  shelf,  causing  sheltered  conditions 
near  the  coast.  In  this  Chapter,  a  preliminary  analysis  of  the  observed  variations  in  Hs, 
0meaa ,  and  <7$  across  the  shelf  is  presented. 

A.  EASTERLY  SWELLS 

The  cross-shelf  variations  for  the  six  easterly  swell  cases  (Table  2)  are  displayed 
in  Figures  9  and  10.  In  all  six  cases,  a  large  decrease  in  wave  heights  across  the  shelf  is 
observed,  with  Hs  at  buoy  XI  on  the  inner  shelf  that  are  between  33  and  53%  smaller 
than  Hs  values  measured  with  buoy  X6  on  the  outer  shelf.  These  decreases  in  wave 
height  are  much  larger  than  the  variations  resulting  from  shoaling  and  refraction  (Herbers 
et  al.,  2000),  suggesting  that  the  strong  attenuation  of  swell  is  caused  by  bottom  friction. 
The  strongest  decay  is  observed  for  the  most  energetic  swells  generated  by  Hurricanes 
Floyd  (September  14)  and  Gert  (September  19  and  20).  Although  the  overall  decay  is 
comparable  for  these  two  events,  the  cross-shelf  evolution  is  different.  During  Hurricane 
Floyd,  strong  initial  swell  decay  (~37%  decrease  in  Hs)  is  observed  across  the  outer  shelf 
between  buoys  X6  and  X5,  followed  by  weaker  variations  across  the  mid-  and  inner 
shelf.  During  Hurricane  Gert,  a  dramatic  decrease  in  Hs  (between  25  and  37%)  is 
observed  on  the  inner  shelf  between  buoys  X3  and  X2,  with  weaker  variations  offshore  of 
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buoy  X3  and  inshore  of  buoy  X2.  The  remaining  three  cases  show  a  more  gradual  wave 
decay  across  the  shelf. 

The  easterly  swell  cases  are  characterized  by  bulk  6mtm  ranging  from  050°  to  135° 
T  (Figure  10a).  The  coastline  and  continental  shelf  break  are  oriented  at  ~  070°  T  and 
100°  T,  respectively,  and  thus,  the  observed  #mean  correspond  to  small  (<  45°)  angles  of 
incidence.  Variations  of  0mean  across  the  shelf  are  small,  except  for  the  more  southerly 
Floyd  and  Gert  cases  that  show  a  ~  30°  shift  toward  normal  incidence.  In  most  cases,  the 
bulk  erg  remained  either  fairly  steady  (between  35°  and  28°)  or  show  a  slight  narrowing 
across  the  shelf  that  is  likely  due  to  weak  refraction  effects.  The  cases  associated  with 
Hurricane  Gert  with  narrow  offshore  spectra  ( ag  «  20  °)  show  a  sharp  increase  of  oe  of 
30%  across  the  inner  shelf  (buoy  X3  to  XI)  (Figure  10b),  possibly  due  to  wave  scattering 
from  bottom  irregularities. 

B.  LARGE  OBLIQUE  ANGLE  SWELLS 

In  the  two  oblique  angle  swell  cases  shown  in  Figure  11,  swells  approach  the 
shelf  break  at  near  grazing  (»  65°  relative  to  normal  incidence)  northeasterly  (December 
11)  and  southeasterly  (October  23)  angles.  In  both  cases  the  significant  wave  height,  Hs 
decreased  by  about  50%  across  the  shelf  (i.e.,  energy  reduction  of  a  factor  of  4)  (Figure 
1  la).  This  decay  is  likely  caused  by  a  combination  of  refraction  and  bottom  friction. 

In  both  oblique  angle  cases,  refraction  effects  are  evident  in  the  gradual  turning  of 

<9mean  across  the  continental  shelf  toward  smaller  oblique  angles  (Figure  1  lb).  However, 
the  observed  refraction  of  the  southeasterly  swell  from  Hurricane  Jose  is  much  stronger 
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(from  165°  T  at  buoy  X6  to  105°  T  at  XI)  than  in  the  northeasterly  swell  case  where  the 
change  is  less  significant  with  swells  arriving  (from  035° T  at  buoy  X6  to  055°T  at  XI). 
These  direction  shifts  are  qualitatively  consistent  with  refraction  computations  presented 
in  Herbers  et  al.  (2000,  cf.  their  Figure  3d). 

The  <re  (Figure  11c)  observed  in  the  northeasterly  case  is  fairly  uniform  across  the 
shelf  (about  30°  to  40°).  In  contrast,  the  southeasterly  swells  are  broader  at  the  offshore 
buoys,  X5  and  X6(aa  «  50-55°)  narrowing  to  38°  at  the  inner  shelf  buoy  XI .  The 

unusually  large  <re  offshore  does  not  appear  to  be  related  to  the  shelf  topography  (buoy 
X6  is  in  deep  water)  and  suggest  possible  scattering  effects  of  the  Gulf  Stream. 

C.  STORM  WAVES 

The  cross-shelf  evolution  of  seas  generated  by  Hurricanes  Floyd,  Irene  and  the 
local  November  2  storm  are  presented  together  in  Figure  12.  The  dominant  waves 
arrived  from  south-southeast  directions  in  all  three  cases,  indicating  fetch-limited 
conditions  on  the  shelf.  The  nor’easter  storm  is  an  unlimited  fetch  case  with  waves 
arriving  from  unobstructed  north-northeast  directions.  Therefore,  the  nor’easter  case  is 
presented  separately  in  Figure  13. 

In  the  three  fetch-limited  cases,  the  observed  wave  heights  decrease  across  the 
continental  shelf  by  about  40-50%  between  buoys  X6  and  XI  (Figure  12a).  These  Hs 
reductions  on  the  inner  shelf  are  comparable  to  the  decay  observed  in  swell  conditions 
(cf.  Figures  10, 12a),  but  the  mechanisms  may  be  quite  different  as  waves  are  actively 
generated  on  the  shelf.  In  all  three  cases,  the  0mean  estimates  show  waves  arriving  from 
south-southeasterly  directions  at  the  outer  shelf  buoys,  X6  and  X5,  shifting  by  about  60- 
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65°  to  near  shore-normal  directions  at  buoy  XI  (1 15-120°  T).  The  large  differences 
observed  between  the  dominant  wave  and  wind  directions  on  the  inner  shelf  (e.g., 
southerly  winds  and  easterly  waves  during  Floyd)  can  be  explained  by  the  extremely 
limited  fetch  for  southerly  winds  and  the  strong  refraction  of  waves  propagating  nearly 
parallel  to  the  local  depth  contours  (Figure  12b).  In  all  three  cases  ae  «  25-30°  at  the 
outer  shelf  buoys  (X4-X6),  but  the  estimates  diverge  on  the  inner  shelf  (Figure  12c).  In 
Hurricane  Irene  ad  increases  to  45°  at  buoy  XI,  whereas  the  cre  values  of  Hurricane 
Floyd  and  the  local  November  storm  case  decrease  slightly  to  20-25°  at  buoy  XI .  These 
differences  are  discussed  in  Chapter  VI.  Variations  of  Hs,  dmt3D ,  and  ae  observed  in  the 
nor’easter  case  are  small  (Figure  13).  Hs  decreased  by  only  25%  across  the  continental 
shelf.  The  observed  0mam  are  about  020°-025°  T  at  buoys  X6  and  X5,  shifting  to  040°- 
045°  T  at  buoys  X4  to  XL  The  u6  are  nearly  uniform  across  the  buoy  transect  varying 
from  25°-30°.  These  weak  variations  indicate  that  refraction  is  of  lesser  importance  in  a 
nor’easter,  and  that  the  wave  field  is  fully  developed  in  equilibrium  with  the  wind. 


26 


VI.  FREQUENCY  DIRECTIONAL  SPECTRA 


The  evolution  of  the  frequency-directional  wave  spectrum,  S(0;  f) ,  across  the 
continental  shelf  is  examined  in  detail  for  the  four  storm  wave  events:  Hurricanes  Floyd, 
Irene,  the  November  2  storm,  and  the  nor’easter. 

Estimates  of  E(f),  0meaa  (/) ,  and  cr0(f )  at  the  peak  of  Hurricane  Floyd  are 
shown  in  Figure  14.  The  frequency  spectra,  E(f),  (Figure  14a)  show  a  large  decrease  in 
spectral  levels  between  offshore  buoy  X6  and  inner  shelf  buoy  XI  in  the  frequency  range 
0.08  to  0.12  Hz  that  corresponds  to  the  peak  of  the  offshore  spectrum.  The  spectral  peak 
level  decreases  monotonically  from  the  outer  shelf  to  the  inner  shelf  (about  a  factor  of  5 
between  buoys  X6  and  XI),  even  though  the  maximum  winds  of  35  m/s  were  observed  at 
the  FRF  station  close  to  the  inner  shelf  buoys  (Figure  5a).  At  both  higher  frequencies  and 
lower  frequencies,  the  spectrum  is  nearly  uniform  across  the  shelf.  The  decrease  in  wind 
sea  peak  energy  on  the  inner  shelf  may  be  caused  by  bottom  damping  and/or  limited 
fetch. 

The  #mean  (/)  estimates  show  large  cross-shelf  variations  caused  by  refraction 
(Figure  14b).  At  the  outer  shelf  buoys,  X5  and  X6,  southerly  angles  (180°  T  through 
200°  T)  are  observed  at  all  frequencies.  Between  buoys  X5  and  XI,  a  shift  in  direction  is 
observed  that  is  largest  for  the  lower  frequencies.  This  is  qualitatively  consistent  with 
refraction  that  is  strongest  for  low  frequency  waves  that  feel  the  bottom  over  the  entire 
shelf.  For  example,  at  buoy  XI,  the  mean  direction  at  high  frequencies  is  about  190°  T 
(close  to  the  local  wind  direction),  and  decreases  with  decreasing  frequency  to  120°  T 
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(near  normal  incidence)  at  the  spectral  peak.  Topographic  effects  play  an  important  role 
in  the  evolution  of  lower  frequency  waves,  but  have  little  effect  on  the  high  frequency 
waves  that  are  in  equilibrium  (i.e.,  saturated)  with  the  local  wind. 

For  all  six  buoys,  crg(f )  is  narrowest  at  the  spectral  peak  (15°-25°),  and  broadens 
at  both  higher  and  lower  frequencies  (Figure  14c).  No  systematic  variations  of  cr0(f ) 
across  the  shelf  are  observed,  but  at  most  frequencies,  the  spreads  are  generally  narrowest 
at  the  inner  shelf  buoys  (X1-X3)  and  widest  at  the  outer  shelf  buoys  (X4-X6). 

Although  Hurricane  Irene  passed  to  the  southeast  of  the  buoy  transect,  producing 
different  conditions  on  the  shelf  compared  to  Hurricane  Floyd,  there  are  many  similarities 
in  the  spectra.  The  spectral  levels  were  not  as  high  as  during  Floyd,  but  the  variations  of 
spectral  levels  across  the  shelf  are  similar  with  about  a  factor  of  5  reduction  between  X6 
and  XI  at  the  spectral  peak  and  weak  variations  at  lower  and  higher  frequencies. 

At  low  frequencies,  the  0mean  (/)  estimates  of  Irene  are  similar  to  those  observed 
during  Floyd,  showing  angles  of  refraction  from  about  180°  to  120°  at  inner  shelf  buoy 
XI  (Figure  15b).  At  high  frequencies  (>  0.15  Hz)  a  wide  range  of  directions  is  observed. 
At  XI,  high-frequency  waves  arrive  from  the  northeast  while  at  X6,  south-southeasterly 
angles  are  observed  (Figure  15b).  These  differences  can  be  explained  based  on  the  local 
wind  observations.  At  the  time  of  the  analysis,  Irene  was  just  passing  to  the  south  of  the 
buoy  transect  and  the  winds  were  rapidly  shifting  from  south  to  north  at  the  FRF  station, 
while  wind  measured  at  the  NOAA  buoy  remained  in  the  southerly  direction  until  much 
later  in  the  day  on  October  1 8  (Figure  5b).  The  high  frequency  part  of  the  spectrum 
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apparently  responds  rapidly  to  the  sudden  shift  in  wind  direction  on  the  inner  shelf, 
resulting  in  a  shift  to  northeasterly  directions  at  buoys  XI  and  X2. 

The  corresponding  cre(f)  estimates  show  a  large  increase  in  spreading  at  high 
frequencies  between  the  outer  shelf  (30-45°)  and  the  inner  shelf  (55-70°)  that  is  likely 
caused  by  bimodal  directional  spectra  on  the  inner  shelf  with  new  northeasterly  seas  and 
old  southeasterly  swells.  At  lower  frequencies,  the  &#(/)  estimates  are  similar  to  those 
observed  in  Floyd  with  smaller  values  at  the  spectral  peak  frequency  and  weak  variations 
across  the  shelf. 

In  the  November  2  storm  case,  the  energy  spectrum  decreases  uniformly  across 
the  shelf  (Figure  16a).  At  all  frequencies,  spectral  levels  at  buoy  XI  are  almost  an  order 
of  magnitude  smaller  than  those  observed  at  X6.  This  large  decrease  is  consistent  with 
large  variations  of  winds  across  the  shelf.  The  outer  shelf  experienced  much  stronger 
southerly  winds  (19  m/s  at  the  NOAA  Buoy)  than  the  inner  shelf  (9  m/s  at  the  FRF) 
(Figure  6a).  Estimates  of  dmam  (/)  (Figure  16b)  show  a  large  shift  from  southwesterly 
directions  on  the  outer  shelf  to  southeasterly  directions  on  the  inner  shelf.  At  buoy  X6, 
high  frequency  components  propagate  in  roughly  the  local  southerly  wind  direction.  But 
at  buoy  XI,  the  local  wind  (200°  T,  Figure  6a)  and  high  frequency  wave  directions 
(135°T,  Figure  16b)  differ  by  75°.  These  observations  indicate  that  high  frequency 
waves  observed  on  the  inner  shelf  are  generated  by  stronger  winds  on  the  outer  shelf  and 
refracted  over  the  sloping  shelf  bathymetry.  At  high  frequencies,  the  cre(f)  estimates 
show  a  systematic  decrease  across  the  shelf  (Figure  16  c)  with  much  narrower  spread  at 
buoy  XI  (25°-30°) ,  compared  with  buoy  X6  (45°-50°),  that  is  qualitatively  consistent 
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with  the  expected  directional  narrowing  of  waves  refracting  from  grazing  angles  toward 
normal  incidence. 

The  nor’easter  case  on  November  1 1  shows  weaker  variation  across  the  shelf 
compared  to  the  previous  cases.  At  high  frequencies,  spectral  levels  are  uniform  across 
the  shelf,  characteristic  of  a  saturated  spectral  tail  in  equilibrium  with  the  wind.  At  low 
frequencies  (<0.15  Hz)  spectral  levels  decrease  monotonically  across  the  shelf, 
suggesting  attenuation  of  long  waves  by  bottom  friction  (Figure  17a).  At  the  peak 
frequency  (0.1 1  Hz),  the  spectral  decay  between  buoys  X6  and  XI  is  about  a  factor  of  2. 
At  all  buoys,  <9mean  (/)  is  within  the  range  15°-45°,  close  to  the  wind  direction,  except  at 
very  low  frequencies  (cf.  Figures  17b  and  6b).  The  cre(f)  estimate  is  very  similar  at  all 
six  buoys,  with  smallest  values  of 20-25°  at  the  spectral  peak,  broadening  at  lower  and 
higher  frequencies  (Figure  17c). 

The  directional  spread  observed  for  all  twelve  cases  is  summarized  in  Figure  18. 
The  common  trend  for  all  buoys  is  a  minimum  spread  at  the  spectral  peak  (/ /  fp=  1),  and 
increasing  spread  for  both  lower  (f!fp<\)  and  higher  (f/fp  >1)  frequencies.  These 
observations  are  consistent  with  deep  water  observations  by  Mitsuyasu  et  al.  (1975)  and 
Hasselmann  et  al.  (1980).  The  range  of  spread  values  at  f=fp  (listed  in  Table  3)  does  not 
vary  much  between  the  outer  and  inner  shelf  buoys.  At  higher  frequencies  the  observed 
directional  spreads  are  more  variable  on  the  inner  shelf  (cf.  buoys  XI  and  X6),  reflecting 
the  complex  wind  generation  (i.e.,  limited  fetch)  and  topographic  effects  (e.g.,  strong 
refraction  along  coast  hugging  waves,  scattering  effects)  near  the  coast  that  can  both 
increase  or  decrease  the  directional  spread. 


30 


Buoy 

Directional  Spread  Range  at  fp 

XI 

—  13°  to  34° 

X2 

14°  to  28° 

X3  ■  1 

12°  to  29°  i 

X4  ’ 

~~  10°  to  28° 

X5 

12°  to  27° 

X6 

12°  to  29° 

Table  3.  Directional  spread  range  for  buoys  XI  to  X6  at  f/fp=  1  in  Figure  18. 
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rv.  CONCLUSIONS 


During  the  SHOWEX  (SHOaling  Waves  Experiment),  a  cross-shelf  transect  of 
six  Datawell  Directional  Waverider  buoys  was  deployed  on  the  North  Carolina 
continental  shelf  to  investigate  the  transformation  of  wind-sea  and  swell  spectra  in 
shallow  water.  A  wide  range  of  wave  conditions  was  observed  during  the  three  month 
deployment  including  several  major  hurricanes.  A  preliminary  analysis  of  the  buoy 
observations  is  presented  here.  Estimates  of  the  significant  wave  height,  Hs,  frequency 
spectrum,  E(f),  mean  direction,  0mean(/)5  and  directional  spread,  c»e(/)  were  extracted 
from  the.buoy  measurements  using  standard  techniques.  The  evolution  of  wave 
properties  across  the  shelf  is  examined  to  evaluate  the  relative  importance  of  bottom 
friction,  refraction  effects,  and  limited  wind  fetch. 

The  analysis  is  focused  on  twelve  case  studies  that  include  remotely  generated 
swell  from  Hurricanes  Floyd,  Gert  and  Jose,  energetic  wind  sea  conditions  during  the 
passage  of  Floyd  and  Irene,  and  a  typical  nor’easter  for  this  region.  In  all  swell  cases,  a 
large  (33-53  %)  decrease  in  wave  height  across  the  shelf  was  observed  that  confirms  the 
dominant  role  of  bottom  friction  in  the  transformation  of  swell  across  a  wide  continental 
shelf.  These  results  are  consistent  with  earlier  observations  at  the  same  site  (Herbers  et 
al.,  2000). 

Estimates  of  0mean(/)  show  the  gradual  refraction  of  waves  toward  shore-normal 
directions.  For  swells  arriving  at  the  shelf  break  at  a  large  oblique  southerly  angle,  0mean 
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shifts  of  up  to  90°  are  observed  on  the  inner  shelf,  enhancing  the  wave  height  decay 
across  the  shelf. 

During  the  passage  of  Hurricane  Floyd,  wave  heights  decreased  from  7.8  m  at  the 
shelf  break  to  4.9  m  on  the  inner  shelf  even  though  the  strongest  winds  were  observed 
near  the  coast.  Possible  mechanisms  that  contribute  to  the  wave  height  decrease  are 
enhanced  bottom  friction  in  shallow  water,  limited  fetch  on  the  inner  shelf  for  southerly 
winds  and  strong  refraction.  Similar  evolution  was  observed  during  Hurricane  Irene.  In 
contrast,  weak  wave  heights  decay  (23%  decay)  between  the  offshore  and  inner  shelf 
buoys)  was  observed  during  a  nor’easter  storm.  Waves  are  more  uniform  in  this  case 
because  the  fetch  is  unlimited  for  northeasterly  winds  and  refraction  effects  are  weaker 
for  waves  arriving  from  northeasterly  directions. 

For  all  twelve  cases  at  all  six  buoys,  the  ce(/)  estimates  consistently  show  a 
minimum  (25°  to  30°)  near  the  spectral  peak  frequency.  Cross-shelf  variations  of  gq  are 
generally  small,  with  the  exception  of  the  directionally  narrowest.  Hurricane  Gert  swell 
observations  that  show  a  dramatic  broadening  across  the  inner  shelf  possibly  owing  to 
bottom  scattering  effects. 
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APPENDIX.  FIGURES 
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Latitude 


Figure  3.  Diagram  of  a  Datawell  Directional  Waverider  buoy  and  mooring  (Datawell  bv, 
1995). 
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Figure  4.  Hurricane  tracks  for  Floyd,  Gert,  Irene  and  Jose  in  1999. 
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Figure  5.  Time  series  of  wind  and  wave  observations  during  a)  Hurricane  Floyd: 
September  15-18, 1999,  and  b)  Hurricane  Irene:  October  17-19, 1999.  Wind  speed 
observations  from  the  FRF  station  and  NOAA  buoy  44014  (top),  wind  direction  (middle) 
and  significant  wave  height,  Hs  at  the  offshore  buoy,  X6  (bottom).  The  vertical  lines 
indicate  the  time  period  analyzed  for  each  event. 
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Figure  6.  Time  series  of  wind  and  wave  observations  during  a)  a  storm  on  November  2- 
4, 1999,  and  b)  a  nor’easter  storm  on  November  10-13, 1999.  Wind  speed  observations 
from  the  FRF  station  and  NOAA  Buoy  44014  (top),  wind  direction  (middle),  and 
significant  wave  height,  Hs,  at  the  offshore  buoy,  X6  (bottom).  The  vertical  lines  indicate 
the  time  period  analyzed  for  each  event. 
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a 


Bulk  Theta  Mean  (al ,  bl ) 

b 


Bulk  Directional  Spread  (al,  bl) 

Figure  8.  Scatter  plots  of  the  second-order  (a2,  b2)  vs  first-order  (al,  bl)  approximations 
of  a)  bulk  mean  propagation  direction,  0mean,  and  b)  bulk  directional  spread,  ctq. 

Estimates  are  compared  for  the  12  case  studies  listed  in  Table  2. 


43 
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Ge  vs  longitude  for  large  oblique  angle  swell  cases. 
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Bulk  Directional  Spread  (deg)  Bulk  Mean  Direction  (deg) 


Figure  13.  a)  Hs,  b)  0mean,  and  c)  ce  vs  longitude  for  nor’easter  swell  case. 
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Figure  14.  Frequency  and  direction  spectra  of  all  6  buoys  (XI -X6)  for  Hurricane  Floyd, 
a)  The  energy  density  spectra,  E(f);  b)  mean  propagation  direction,  0mean  (/);  and 
c)  the  directional  spread,  Ge(/). 
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J?  tene:  1 0/1 8  00  EST  to  1 0/1 8  03  EST 


Figure  15.  Frequency  and  direction  spectra  of  all  6  buoys  (XI -X6)  for  Hurricane  Irene, 
a)  The  energy  density  spectra,  E(f);  b)  mean  propagation  direction,  Qmean  (/);  and  c)  the 
directional  spread,  ae  (/). 
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Directional  Spread  (deg)  Directional  Spread  (deg)  Directional  Spread  (deg) 


Figure  18.  ere  estimates  of  all  twelve  case  studies  vs  normalized  frequency  (f/fp )  for  each 
buoy. 
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